Long ranged singlet proximity effect in ferromagnetic nanowires 
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Recently a long ranged superconductor/ferromagnet (S/F) proximity effect has been reported in 
Co crystalline nanowires 0, Nature, 6 389 (2010)]. Since the authors of [l[ take care to avoid the 
existence of magnetic domains, the triplet character of the long ranged proximity effect is improbable. 
Here we demonstrate that in the one-dimensional ballistic regime the standard singlet S/F proximity 
effect becomes long ranged. We provide an exact solution for the decay of the superconducting 
correlations near critical temperature (T c ) and for arbitrary impurities concentration. In particular 
we find a specific regime, between the diffusive and ballistic ones, where the decay length is simply 
the electronic mean-free path. Finally possible experiments which could permit to elucidate the 
nature of the observed long ranged proximity effect in Co nanowires are discussed. 

PACS numbers: 74.45. +c, 74.78. Fk 85.25.Cp. 



It is well known that the superconducting proximity ef- 
fect in a diffusive ferromagnetic (F) metal is rather short- 
ranged due to the large value of the ferromagnetic ex- 
change field E ex ~ (500 — 5000) K acting on the electron 
spins. In the fcrromagnet, the singlet Cooper pair wave 
function experiences an oscillatory decay whose charac- 
teristic length £f = \J HD / E ex enters the nanoscopic 
range for typical values of the diffusion constant D and 
exchange field in the ferromagnet [H, Q. This is in a 
sharp contrast with the corresponding decaying length 
in normal metals = y^hD/ksT ^ £,f which can reach 
microns at low enough temperatures. 

The presence of the non-uniform magnetization may 
strongly modify the S/F proximity effect for singlet su- 
perconductors, and in the diffusive regime the induced 
triplet correlations can penetrate at the large distance £„ 
In the clean limit and in the presence of domains, the 
triplet proximity effect also increases the decay length up 
to the thermal decay length £ n& = hvp/ksT where 
vf is the Fermi velocity. The long range proximity ef- 
fect has been observed in Ho Q and in the half metal 
Cr02 0- Probably its origin is related with the presence 
of magnetic domains in Ho, and with magnetic scatter- 
ing at the S/F interface in the case of Cr02. Recently 
the controlled transition between the short ranged singlet 
and the long ranged triplet proximity effect has been re- 
orted in Co based Josephson junctions using Co/Ru/Co 
and Ho/Co/Ho Q as magnetic barriers, respectively 
Previously the 200 nm decay length for the coherence ef- 
fect was observed in a system with mesoscopic Co wire 
lOf. The geometry of the Co electrodes in the setup of 
10] was nevertheless 2D-like and probably the observed 
long ranged proximity effect was related with the pres- 
ence of the magnetic domains which generated the triplet 
superconducting correlations 



Very recently the Josephson current through the Co 
crystalline nanowires as long as 600 nm was reported 
[if. Interestingly the experiments on the Co crystalline 
nanowires [l[ were performed after applying a strong 
magnetic field which makes improbable the existence 
of the magnetic domains. The diameters of the Co 
nanowires |l| were 40 nm and 80 nm and the proxim- 
ity effect was substantially weaker for thicker nanowires. 
Note that in the S/F/S structures with Co film as a 
F layer, the long ranged proximity effect is absent and 
the characteristic decay length of superconductivity is 
around 5 nm 

mm. Theref ° re the 

nanoscopic charac- 
ter of the Co wire seems to be crucial for the observation 
of the long ranged proximity effect. 

In this Rapid Communication we present a theoretical 
analysis of the S/F proximity effect in the ID case tak- 
ing in mind a thin Co nanowire. For this simple model 
it is possible to obtain an exact solution of Eilenberger 
equations with arbitrary impurity scattering. The decay 
of the superconducting correlations is governed by the 
electron mean free path I. We believe that for the Co 
nanowire I may be in the range (50 — 100) nm which 
could explain the long ranged proximity effect observed 
in [H . We also discuss the alternative scenarios including 
dead layers or triplet correlations. Finally we suggest fu- 
ture experiments to discriminate between those different 
types of explanations. 

The diffusive regime for the proximity effect in ferro- 
magnets corresponds to the limit E ex r <C 1 [l3| which 
could be realized in weak ferromagnets with CuNi or 
PdNi alloys as F layer mi lh but not in the systems 
with strong ferromagnet like Co, Ni or Fe 111 Il2j. How- 
ever in ballistic regime the S/F proximity effect is sub- 
stantially different. The Cooper pair wave function oscil- 
lates inside the ferromagnet at the characteristic length 
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fb 



HvF/E ex [ly] but it decays much slowly : in 3D 
system as ~ (£f b /x) , where x is the distance from the 

S/F interface, in 2D systems as 
decay in ID case [17 



(€fb/ x ), and no 
1 1 Si ] - Sure this statement is true 
for the distances x smaller than thermal decay length 
£ rl6 = HvF/ksT and mean free path I = vft, where r 
is the average time between impurity scattering events. 
Otherwise the range of the proximity effect will be fixed 
by the smallest length 



Near the superconducting critical temperature, the 
Eilenberger equations [22j may be linearized and they 
read for the ID case: 
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where /+ and /_ are the anomalous Eilenberger Green 
functions for vp > and vf < respectively and w is the 
fermionic Matsubara frequency. Note that (/+ + /-)/2 
is simply the average over the Fermi surface in the ID 
case. 

After introducing the symmetric and antisymmetric 
combinations F s = (/+ + /_)/2 and F a = (/+ - /_)/2, 
one obtains that F s obeys the following second order dif- 
ferential equation: 



(cu + iE ex )F s 



d 2 F s 



iE ex + — 
It 



dx 2 



while the antisymmetric combination can be obtained 
from 

F. = -"i (. + . E „ + ^'§. (3, 

Let us now sketch two limiting cases of Eq.Q. In the 
diffusive limit E ex r <C H, Eq.Q becomes: 



(w + iE ex )F s 



vpr d 2 F s 
2 dx 2 







(4) 



which is the Usadel equation with the diffusion coeffi- 
cient D = vpr. In the ballistic limit E ex r 3> ft, Eq.Q 
becomes: 



(uj + iE ex ) F s 



d 2 F s 



4 dx 2 



(5) 



which contains the same information as the Eilenberger 
equation Eq.((T]) in the absence of disorder (1/r = 0). 

Solving the two previous equations (|4I5|) for E ex ^> T c 
shows that in the diffusive regime the characteristic os- 
cillatory and damped behavior of F s occur with the same 
characteristic length £ ^ , whereas the ballistic regime may 
exhibit no attenuation at all. In the following, we inves- 
tigate the effect of a small amount of impurities into the 



ferromagnetic wire. We show that this regime is char- 
acterized by two distinct lengths: the ballistic ferromag- 
netic oscillation one £f b and the damped one I induced 
by the impurities into the wire. 

The general solution of Eq.Q for arbitrary impurity 
scattering is simply: 



F s = Ae qx + Be- qx , 



(6) 



with 
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1" = — 



1 uj + iE ex + j- ) (u + iE ex ) . (7) 



Using Eq.Q and the approximation E ex r 3> 1, one ob- 
tains the antisymmetric anomalous propagator: 



F a 



-Ae qx + Be- gx . 



(8) 



Assuming superconducting electrodes much thicker than 
the F nanowire, we apply the rigid boundary condition 
f + (x = -L/2) = Ae-^^/u). Then we find the Eilen- 
berger propagator for the superconducting correlations 
originating from the left electrode: 



/-+ 



\ e -i V /2 e -q(x+L/2) 



(9) 



where we suppose the temperature close to T c . The sim- 
0, (2) ilar condition /_ (a; = L/2) = Ae"^/ 2 /w yields the cor- 
relations from the right superconductor: 



f_ = _ e iv/2 e «(x-L/2)^ 

UJ 



(10) 



We now evaluate the Josephson supcrcurrent for this sin- 
gle channel situation: 



v^/A\ 2 (2E^i 
= 2nei>oVFT I — I cos 



VF 
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(11) 



where isq = kpl^jEp = (Efcl)~ x is the one-dimensional 
density of states per unit length (including 2 spin direc- 
tions). 

In the pure limit, I 3> L, the critical current exhibits 
undamped oscillations: 



v^/A\ (2E ex L 
I c o = 2tt ev qv pT \ — ) cos 
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(12) 



whose amplitude gives an estimate for the single channel 
critical current. The above result was derived near the 
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critical temperature. Concerning order of magnitudes, it 
can be safely extrapolated to low temperature by substi- 
tuting A = A(T) by A(T = 0) ~ T c . We introduce the 
quantum of resistance Rn = h/e 2 and find: 

IcO ~ ^eiy v F — ~ (13) 



ei? 



Hence we can estimate for the experimental conditions 



1]: 



cRnI c q ~ 1 mV and so J c , 



40 nA. 



(14) 



This maximal current per mode is damped by the expo- 
nential factor e~ L l 21 . 

For a metallic nanowire with a cross section S, the 
number M of transverse channels typically amounts 
S/a 2 ~ (d/a) 2 ~ (10 4 - 10 5 ). This relatively large num- 
ber of channels may compensate the exponential sup- 
pression of the critical current which can be evaluated 
as I c = Me~ L l 21 I c q. In experiment [![, the critical cur- 
rent exceeds 10 fiA and M ~ 10 5 (d = (40 - 80) nm, 
L = 600 nm and a = 0.1 nm) thereby requiring a rather 
long mean free path, namely I > 60 nm. On the basis of 
the Drude model [l[ , the estimate of the mean free path I 
is quite sample dependent with values ranging in the win- 
dow I ~ (2 — 10) nm. In our approach such short mean 
free paths would lead to a very strong suppression of the 
current. Nevertheless we might question if those values 
I ~ (2 — 10) nm represent the genuine mean free path 
inside the ferromagnet. Indeed the authors [lj note an 
important diffusion of the W atoms inside Co nanowires. 
Therefore the region near the W/Co interface may con- 
tribute to the measured resistivity and the genuine or 
intrinsic mean-free path I (defined between the contact 
regions) of the Co nanowire may be substantially higher 
than the simple Drude estimate. 

Long range triplet proximity effect, caused by some 
magnetic inhomogencity, is also a possible scenario to 
explain the experiment |l|. A novel kind of Joscphson 
junctions S/Xi/F/X 2 /S with 3 magnetic layers (Xi, F, 
X2) has been proposed la] and successfully implemented 
in a controlled way d, Q . A central region (F) is used 
to suppress singlet Josephson current while the side do- 
mains create triplet correlations. In the original theoret- 
ical set-up, the Xi, F and X2 were three monodomains 
with noncollinear magnetizations. In the recent exper- 
iment [U, Xi and X2 were layers of weak ferromagnets 
while F was itself a Co bilaycr, Co/Ru/Co. In Ref.[9( the 
synthetic Ho/Co/Ho ferromagnetic trilayer was used to 
the same purpose. More specifically it was predicted [f| 
and confirmed experimentally Q that the production of 
triplet correlations is maximal for an optimal size of the 
Xi and Xo layers, which is on the order of £j h . The au- 
thors of [1| mention that though the Co nanowires must 
be single domain, the inhomogeneous magnetic moments 
may be produced in the W-Co contact regions. Follow- 
ing [j| and [H, the size of such regions must be of 



order to maximize the triplet long range effect. Nev- 
ertheless this contact region is likely to be of the atomic 
size a and in such a case the triplet Josephson current 
would be reduced by a factor (a/^ b ) 2 ~ (l0~ 2 — 10 _1 ). 
The optimal triplet supercurrent is M ~ 10 5 times the 
single channel S/N/S supercurrent (0.04 fj,A). Hence our 
estimation of the triplet effect in Q is I c = (4 - 40) nA 
which is still in agreement with the reported supercur- 
rents. 

In order to discriminate between the two scenarios 
described above, it is crucial to know whether the Co 
nanowires are in the diffusive or in the ballistic limit, 
and also to have an experimental determination of the 
Josephson current-phase relation 0, Q ■ Indeed the pres- 
ence or absence of a second harmonic at low temperature 
allows a distinction between those two scenarios. In the 
singlet correlation scenario, the second harmonic will be 
strongly suppressed by disorder (by an additional factor 
e -L/2i re i a ti ve to the first harmonic), while for triplet 
proximity effect the second harmonic is expected to be 
of the same order of magnitude as the first one, at low 
temperature. 

Besides, it would also be interesting to know if the 
S/F/S structures with Co nanowire may be in the it state. 
This can be determined by realizing a SQUID with the Co 
nanowire Josephson junction as one arm and a standard 
Josephson junction on the other arm. 

Note that a third possible scenario might involve a 
dead layer at the surface of the Co nanowire, where mag- 
netism is strongly suppresses thereby leading to the prox- 
imity effect reported in [l]]. Indeed the Co wires are sur- 
rounded by an insulating oxide coating and a dead layer 
(without ferromagnetic order) might appear between this 
insulating layer and the Co wire core. The presence of 
such dead layer was already signaled in Co films, with 
typical size 1 nm [ll[. In the extreme case, the mag- 
netism is completely absent within this layer and the cor- 
responding supercurrent would not have the exponential 
decay. In such a case the number of channels would be 
Mgurf = d/a ~ 10 3 and the critical current would reach 
quite large values 10 2 fiA, even for long nanowires. Note 
that in presence of such a deadlayer and at low temper- 
ature, the second harmonic of the current-phase relation 
would be as large as the first harmonic, similar to the 
case of the triplet long ranged proximity effect. 

In conclusion, we have found that the long ranged 
S/F proximity effect reported recently in [lj can be ex- 
plained by three distinct scenarios. In the first scenario, 
the Josephson current is associated with weakly damped 
singlet superconducting correlations. In the second sce- 
nario, the contact regions produce triplet correlations 
thereby leading to long ranged proximity effect. Finally 
in a third scenario a non magnetic deadlayer provides 
a channel for long ranged propagation of the supercon- 
ducting correlations. We further suggest to realize the 
low-temperature determination of the current-phase re- 
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lation to discriminate between those scenarios. 
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